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ABSTRACT 


The flutter stability of flap bending, lead-lag bending, and torsion of helicopter 
rotor blades in hover is investigated using a finite element formulation based on 
Hamilton’s principle. The blade is divided into a number of finite elements. Quasi- 
steady strip theory is used to evalua(<;e the aerodynamic leads. The nonlinear 
equations of motion are solved for steady-state blade deflections through an iterative 
procedure. The equations of motion are linearized assuming blade motion to be a 
small perturbation about the steady deflected shape. The normal-mode method 
based on the coupled rotating natural modes is used to reduce the number of 
equations in the flutter eigenanalysis. First the formulation is applied to single- 
load-path blades, for example, articulated and hingeless blades. Numerical results 
show very good agreement with existing results obtained using a modal approach. 
The second part of the application concerns multiple-load-path blades, namely 
bearingless blades. The flexure of a bearingless blade consists of multiple beams 
(flexbeams and torque tube) leading to redundancy. The formulation is modified so 
that the multibeams of the flexure could be modelled individually. Numerical results 
are presented for several analytical models of the bearingless blade. Results are 
also obtained using an equivalent beam appraoch (the common approach) wherein 
a bearingless blade is modelled as a single beam with equivalent properties. The 
comparison between the two sets of results show that the equivalent beam modelling 
is inaccurate. 
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NOMENCLATURE 


a lift-curve slope 

A blade cross-section area 

B\,B2 blade cross-section integrals 

e blade chord 

Cl blade-section lift coefficient 

blade-section pitching moment coefficient 
Ct rotor thrust coefficient 

Ci,C2 blade cross-section integrals 
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cj aerodynamic-center offset from elastic axis, positive aft 
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E Young’s modulus 
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GJ blade elastic torsion stiffness 

Iy,Ig blade cross-section moments of inertia in the flap and lead-lag direc- 

tions, respectively 

polar radius of gyration of blade cross section, \/{Iy + Iz)/A 
km mass radius of gyration of blade cross section, \J{k^^ + 

limit principal mass radii of gyration of blade cross section 
l{ length of ith element 
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blade cross-section air velocity components in negative i| and f direc- 
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induced inflow 
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local axial coordinate of the ith element 
blade cross-section angle of attack 
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virtual work done due to aerodynamic loads 
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respectively 
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rotor blade angular velocity 

aoidz 

d()ldii> 


O' 

0 




1. INTRODUCTION 


Helicopter dynamics is a complex problem involving nonlinear structural, in- 
ertial, and aerodynamic forces. For efficient blade designs, there is an increasing 
trend toward nonuniformity in structural and aerodynamic properties, The root 
geometry is complex. 

Most of the helicopters in service employ either an articulated or a teetering 
rotor. In recent years, there has been considerable interest in alternate designs with 
an eye toward mechanical simplicity and increased maintainability. The hingeless 
and bearingl^s rotors are the outgrowth of such efforts. A conventional articulated 
rotor blade has mechanical flap and lag hinges, a lag damper and a pitch bearing. 
Figure 1 shows the schematic of an articulated blade hub and root. A hingeless 
blade, as the name implies, lacks the hinges and is cantilevered at the hub. It does 
have a pitch bearing for pitch control. A bearingless design eliminates the hinges 
and the pitch bearing as well. Only some hingeless and bearingless designs include 
lag damperb. The bearingless blade has an elastic flexure consisting of flexbeams 
and a torque tube to facilitate pitch changes. An example (Fig. 2) that has been 
tested in flight and wind dunnel is the Boeing* Vertol designed Bearingless Main 
Rotor, BMR (Dixon and Bishop (1980); and Warmbrodt and McCloud (1981)). The 
BMR flexure (Fig. 3) extends from the hub to about 25%blade radius where it 
is connected to the outboard blade by a rigid clevis. The flexure consists of two 
parallel, torsionally-soft, open-section beams and a torsionally-stiff (flexurally-soft) 
torque tube. The inboard end of the torque tube is connected to a pitch link via 
a pitch horn at the hub; and the outboard end is attached to the clevis. The 
pitch control to the blade is applied through the torque tube by rotating the blade- 
root clevis which elastically twists the flexbeams. The blade outboard of clevis 
is of airfoil section and is similar to a conventional blade. The articulated and 
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hingeless blades are referred to as single-load-path blades and the bearingless blade 
as muliiple-load-path blade, This is because the blade loads can be transmitted to 
the hub by redundant paths of the dexure in the case of a bearingless blade. 

To help in the design of helicopter rotor blades, it is necessary to solve for trim 
deflections and to determine the aeroelastic stability in hover and forward flight. 
The aeroelastic stability of articulated blades is well understood, In the case of 
hingeless and bearingless blades, the absence of hinge motion modifies the dynamic 
characteristics of the blade and may significantly influence the aeroelastic stability. 
Not enough is known about the parameters afiTecting the aeroelastic behavior of 
bearingless blades. The influence of multilple load paths on aeroelastic stability is 
not fully understood. 

The simpl^t form of roptor blade representation b the rigid blade model with 
spring restrained hinges [Ormiston and Hedges (1072) and Kaza and Kvaternik 
(1976)]. For this model, the trim solution is generally obtained from coupled linear 
or nonlinear flap-lag equations and the flutter stability is determined assuming linear 
perturbation motion about the trim solution. This type of modelling is satisfac- 
tory for simple blade configurations such as articulated. A more appropriate repre- 
sentation is to treat the rotor blade as an elastic beam. Houbolt and Brooks (1958) 
have derived the differential equations of motion for the combined bending and 
torsion of an elastic blade. Several authors [Hodges and Dowell (1974), Kvaternik 
and Kaza (1976), and Rosen and Friedmann (1978)] have modified these equations 
to consistently include nonlinear structural and inertial terms for moderately large 
deflections. A number of researchers have aualysed the aeroelastic stability of hinge- 
less blades with elastic-beam modelling. Examples are Hodges and Ormiston (1976) 
and Jehnson (1977). The common approach for calculating the trim deflections, 
as well as the aeroelastic staldily, is the modal method (for example, the Galerkin 
method) using coupled natural modes. In general, the trim deflections are assumed 
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to be large and obtained by solving nonlinear steady-state equations; and the flutter 
equations of motion are linearized about the trim state. This linearization simplifies 
the stability analysis. There have been some attempts at determining the stability 
using a complete nonlinear analysis such as the limit cycle analysis. These have 
been for simple blade models. For example, Friedmann and Tong (1973) consider a 
flap-lag blade; and the analysis of Chopra and Dugundji (1979) incl.ides the torsion 
degree of fraeclom but the blade is modelled as rigid with spring restrained hinges. 

The theoretical analysis of a bearit*sless blade is more involved than that of a 
hingeless blade primarily due to the redundancy of the flexure. Several authors have 
examined the aeroelasUc stabililty of bearingless blades. Bielawa (1976) has given a 
theoretical analysis for composite bearingless rotor blades and experimental results 
from dynamically scalec’ models. The formulation is made for a blade with one 
flexbeam and a torque tube. Hodges (1979) has analyzed bearingless configurations 
with single flexbeam. The bending stiffness of the torque tube is neglected. The 
outboard blade is treated as rigid. Hodges has used an equivalent-beam approach 
to model the twinbeam root flexure of the BMR blade. Harvey (1976) has modelled 
a single-flexbeam bearingless blade by lumped massed connected by elastic beam 
elements. If any of the analyses mentiond above are used to study the BIS^ blade, 
then the twinbeams of the the flexure has to be modelled as single flexbeam with 
equivalent properties. The equivalent properties can be obtained by matching the 
fundamental frequencies of the equivalent-beam model with that of the multibeam 
blade. 

With the modal approach, it becomes increasingly difficult to handle geometric 
coomplexities. For example, it is difficult to efectively model the multibeam flexure 
of a bearingless blade. The finite element method [Gallagher (1975)] has been 
extensively used in linear structural analysis. The structure is divided into a number 
of finite elements and the application of energy principles, or the method of weighted 
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residuals, yields approximate expressions for forces (inertial, elastic, etc.) over 
each element, The global equations of motion are obtained by assembling the 
elements. Nonuniform properties can be easily accommodated. The finite element 
method is very flexible and the formulation can be adapted to different rotor blade 
configurations with a few modifications. Multibeams of a bearingless blade can be 
modelled individually. 

The finite element method has been applied for the determiniation of free- 
vibration characteristics of rotating one dimensional (flap bending) beams by many 
authors. A few examples are Nagaraj and Shanthakumar (1975), Putter and Manor 
(1978), and Hodges and Rutkowski (1981). It is common practice to use simple beam 
elements (represented by end nodes) satisfying the continuity of displacement and 
slope between elements. Refined elements with internal nodes may be used instead. 
One can attain a desired accuracy by using either simple or refined beam elements, 
but, if simple elements are used, a greater number of them may be necessary. 
However, the assembly of simple elements results in narrowly banded matrices; the 
assembled matrices for refined elemnts are more populated. Murty and Raman 
(1980) have used the finite element method for computing the nonlinear response 
of a rotating beam under prescribed forces. 

Yasue (1978) has obtained the natural-vibration characteristics of a hingeless 
blade in combined flap bending, lag bending, and torsion using a finite element 
formulation. The normal mode method based on these natural modes is utilized to 
calculate the trim solution and the response of the blade under gust loading. Yasue's 
finite elements exhibit the interelement continuity of the torsional deflection as well 
as its slope; and hence a blade with jumps in torsional rigidity cannot be handled. 
Friedmann and Straub (1980) and Straub and Friedmann and Straub (1981) have 
applied a finite element formulation based on the method of weighted residuals to 
calculate the free-vibration characteristics of a flap-lag hingeless blade. Then, the 
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trim state, flutter stability, and response of the blade in hover and forward flight 
have been computed using a modal approach. These three rdeferences use the finite 
element method just to obtain the natural modes, the trim solution is based on a 
modal approach. 

In the present research, a finite element formulation is developed to investigate 
the flutter stability of flap bending, lead-lag bending, and torsion of rotor blades in 
hover. The formulation is made for a nonuniform blade with pretwist and precone, 
and having chordwise offsets of the center of mass, aerodynamic center, and tension 
center from the elastic axis. The deflections considered are the ivxial, lead-lag, and 
flap deflections of the elastic axis and a torsional deflection about the elastic axis. 
The spatial dependence of the deflections is made discrete by dividing the blade 
into a number of elements. Approximate expressions for element forces are obtained 
applying Hamilton’s principle; and the assembly of elements results in the global 
equations of motion in terms of the nodal degrees of freedom. The aerodynamic 
loads are based on quasi-steady strip theory. The trim solution is computed from 
nonlinear steady state equations without making a modal transformation. The 
coupled rotating natural-vibration characteristics are calculated about the trim 
condition. The flutter motion is assumed to be a small perturbation about the 
trim state. The normal mode method based on the rotating natural modes is used 
to solve the linearized flutter equations as an eigenvalue problem. Flutter stability 
is inferred by studying the flutter eigenvalues. 

Here it should be pointed out that the emphasis of the present study is on the 
feasibility of the application of the finite-element method to rotor blade aeroelas- 
ticity; and hence, importance is not given to certain aerodynamic aspects such as 
compressibility effects and the effects of preceding and returning wakes. A brief 
review of the research in these areas follows. It is popular to use a quasisteady 
strip theory approximation without wake effects for estimating the unsteady loads 
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on a rotor, one of the main reasons being the difficulties involved in accuractely 
representing the Theodorsen function e(/s). The difficulties arise because, unlike 
that of fixed >vings, the wake of one rotor blade interacts with that of the others. 
The first significant attempts at unsteady aerodynamic theories for rotary wings in 
incompressible flow are due to Loewy (1057), Timman and Van De Vooren (1957), 
and Jones (1958). A theory to study the compressiblity effects including wake is 
presented by Jones and Rao (1970); aerodynamic coefficients for certain Mach num- 
bers and frequencies are given in this paper. In another study, on the effects of tip 
vortices on rotor blades in hovering flight, Jones and Rao (1971) colclude that the 
use of two-dimensional strip theory would not result in serious error provided the 
blades oscillate at several cycles per rotation. Murthy and Pierce (1976) formulate a 
technique for the prediction of flutter of a helicopter rotor in hover and axial flight 
in such a manner that various unsteady aerodynamic theories and various types of 
coupling between the degrees of freedom could be accommodated. Kato and Yamane 
calculate the rotor impedence of articulated-rotor helicopters in hover (1979a) and 
forward flight (1979b) using a quasisteady approximation that neglects the effects of 
preceding and returning wakes. In trying to correlate their theoretical results with 
experiments, Kato and Yamane (1981) find that the quasisteady theory does not 
predict the amplitude peaks and clefts which appear near the multiples of the blade 
passage frequency 60 (where h is the number of blades) whereas an analysis with 
Loewy’s function agrees fairly well with the experiment. Various dynamic inflow 
models are included in a flap-lag stability analysis of a rigid blade in forward flight 
by Peters and Gaonkar (1980). Their study concludes that unsteady perturbations 
in the induced flow have a significant effect on both the flap and lead-lag damping, 
particularly for the regressing mode. 

The application of the present method consists of two parts. First, numerical 
results are obtained for uniform single-load-path blades. Results are compared 
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with that of a previous research based on a modal method. Comparisons with 
experimental results are also made. The second part of the application concerns 
multiple-load-path blades. The formulation made for single-load-path blades is 
modified here to include the multibeams of the root flexure of bearingless blades. 
Each of the flexbeams and the torque tube are modelled as individual beams. The 
displacement compatibility conditions, at the clevis, between the inboard beams and 
the outboard beam are used in the modelling of the bearingless blade. Refined finite 
elements with internal nodes are used. Numerical results are obtained for several 
analytical models of the bearingless blade in hover. To determine the accuracy of 
the conventional equivalent- beam approaches, results are obtained for a blade with 
equivalent properties. Comparisons are made between the results of the present and 
equivalent-beam models. 

This is the first attempt to solve the nonlinear trim equations of a flap-lag- 
torsion rotor blade using the finite element analysis directly. This approach enables 
the coupled natural rotaing modes to be calculated about the steady deflected 
positon of the blade. This is the first research to apply the finite element method 
to bearingless blades; and the first one to model the multibiams of the flexure 
individually. 



2. HAMILTON’S PRINCIPLE FORMULATION 

2.1 Coordinate Systems 

The coordinate systems used are shown in Fig. 4. The rotor blade is treated as 
an elastic beam rotating at constant angular velocity n. l^be rectangular coordinate 
system x, y, z is attached to the undeformed blade which is at a precone angle of ^p. 
The origin is at the root of the blade, the z axis coincides with the elastic axis, and 
the y axis is in the plane of rotation pointed towards the leading edge. The deformed 
position of the blade is defined as follows. A point P on the undeformed elastic axis 
undergoes displacements u,v,w in the x,y,z directions, respectively, and occupies 
the position P* on the deformed elastic axis; u is the axial deflection, v the lead-lag 
deflection, and tw the flap deflection. Then the blade cross section containing P* 
undergoes a rotation 0i about the deformed elastic axis. The orientation of the 
deformed-beam cross section with respect to the undeformed- beam cross section is 
described by a sequence of three rotations. This analysis uses the iag-flap-pitch 
sequence of rotations defined in Hodges, Ormiston and Peters (1080). The thrid 
angle in the sequence, is written as 


0 ^ 

and ^ to second order is 


( 2 . 1 ) 



( 2 . 2 ) 


where $ is the pretwist, ^ is the geometric twist, and is the elastic twist due to 
torsion. The orthogonal coordinate system tf, f is attached to the deforemed blade 
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such that the ( axis is tangential to the deflected elastic axis and the r/ and ; axes 
are the principal axes of the cross section. The transformation between the two 
coordinate systems is given os 


/*» 


t-i 

*» 

\*n 

. in 

t 

*» 

1 1 


1 



. *2 , 


(2.3) 


where 7^, and 7*, ly, are the unit vector systems in the two coordinate frames. 
The transformation matrix is defined in Hodges et al (1080) and is reproduced, in 
appendix A. 

2,2 Hamilton’s Principle 

The formulation is based on Hamilton’s principle: 


ft2 

As / {6U-ST-SW)dt^^0 (2.4) 

Jti 

where SU, ST, and 6W are, respectively, the variation of strain energy, the variation 
of kinetic energy, and the virtual work done by external forces. Hodges and Dowell 
(1074) give expressions for SU, and ST. The expression for the variation of the strain 
energy is , 
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ORIGINAL PAGE 13 
OF POOR QUALITY 


5C/- 

+ + EAk\{9 + j 

+ EBx - EB29\(v ^* cos + w” sin 
+ [EC\<i>^' + EC 2 {iv'^ cos9i - a\n9i)]6^*' 

||E7j^v^^costfi + w^Sin^l^— £7A«a^«^ + ~ EB2^^* cos^i 

EIy(^v*^ airx 9 X — ty^costfi^— sintfil^tfv'^ + 

cos ^1 + sin 9x^—EAea^u^ + “ EB2^4>^ *•“ ^1 

- sin 9x — cos cos 9x 


+ 

+ 

+ 


( 2 . 5 ) 


where the axial force F is defined by 


+ k\9^^* — Ca(v^^cos 9x + w^sin tfi)j 


( 2 . 6 ) 


and the section consti</»nts are 
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ORIGINAL PAGE \n 
OF POOR QUALITY 


Am 

ACfi HI 
lym 

/* = 
Ak\ m 

B2 m 

Cim 

Czm 



j l^{v^ + S^)d„di 

/Z 

j j^\-d„d, 
j j^f\d^d, 



( 2 . 7 ) 


The blade cross-section area effective in carrying tension is A. ly and Iz are, 
respectively, the flapwise and chordwise area moments of inertia of the cross section; 
Ak\ is the polar moment of inertia of the cross section; J is the torsional constant 
including cross-section warping; ta is the offset of the tension axis from the elastic 
axis, positive forward; Bi, B2, C\, and C2 are other section constants; X is the warp 
function which is assumed to be antisymmetric about the chord. The cross section 
is assumed to be symmetric about the chord. 

The variation of the kinetic energy for the elastic beam is 
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ST - |m(n2® + 2nu - u)«u 

+ + eg cos tfi) + 2mn(^pti; — ti) 

+ 2me^fl^t)^ cos 9i + V sin j — mi; + meg^ sin 0i 

- m^p(n^x + 2nt) j + mw + mcg^ cos $i 6iv 

“ + fnCl^(^k%^ — k^mi^j cos $i sin $i 

+ meg Q'^ x{‘w* cos 0\ - v' sin 0i ) + meg Q^v sin 0\ 


6v 


+ mCyO^^pX cos — m«j;(i; sin l?l — wcostfi) 

— mcy^n^xcostfi + 2f2tcostfi 

— mcp^n^xsin^l + 2nwsin^l^tfti;^| dx 


S^ 


( 2 . 8 ) 


Where the section integrals are defined by 



( 2 . 9 ) 


The blade mass per unit length is m; eg is the center of mass ofiset from the 
elastic axis, positive forward; and kmi and ^^>2 l>he flapwise and chordwise mass 
moments of inertia per unit length. 

The virtual work 8W of the external forces may be expressed as 
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LySv LmSw + dx 


( 2 . 10 ) 


=* 


where Lu, Lv, Lw, and are the external loads distributed along the length of the 
blade in the axial, lead-lag, flap, and torsion directions respectively and the virtual 
rotation 6 if as defined in Hodges et ol (1680) is 


6if ( 2 . 11 ) 

Since 6U, ST, and as given in Eqs. (2.5), (2.8), and (2.10), are independent of 
the time derivatives of 6u, 6v, Sw, and S^, Eq. (2.4) can be written as 


AH«ir-6r-«Vr-0 (2.12) 

2.3 Elimination of the Axial Displacement 

For single load path blades, the axial displacement u can be eliminated in terms 
of the other deflections v, w, and i> and the centrifugal force F, This is possible 
only for single-load-path blades since the force equilibrium equation in the axial 
direction can be solved a priori. 

Equation (2.12) would result in four equations of motion representing Su, Sv, 
Sw, and 6^, respectively. The Su equation is 

Ftfu' + 2mn0 -I- Lu di *= 0 (2.13) 
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Integrating by parts and noting that 
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■■ 0 at X"mO 


and 


F ••O at X * /?, 
the above equation can be written as 

’i-mn^X’jr2mnv + Lu^Sudx Ml) (2.14) 

Since this equation is true for arbitrary values of Su, 


F^ + mn^x + 2mUv + Lu 0 ' (2.15) 


The external loads of aerodynamic origin and the radial drag force Lu is negih 
gible in hover. So neglecting Lu and integrating Eq. (2.15) results in 



n^x + 2riv^dx 


(2.16) 


The displacement u is eliminated from SU, ST, and SW using Eqs. (2.6) and (2.16). 
2.4 Nondimensionalisation 


After eliminating u, the SU, ST, and SW are nondimensionalised by dividing the 
expressions by where mo is a reference mass per unit length (for example, 
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the value of m at 0.5/?), This is done to see the relative order of magnitude of 
each term. The nondimensional displacements v//?, w/R, and ^ are assumed to 
be of order < where t is a small nondimensional quantity such that < 1. The 
other nondimensional quanities and their assumed orders of magnitude are given in 
Appendix B. The lowest order (i.e, flrst order) terms in 517 and ST are of order c-, 
Terms of order e- and e® are retained while neglecting all terms of order (* (i.e., third 
order) and higher, excepting some linear third order terms which are important for 
the torsion equation. The following simplified nondimensional expressions for SU, 
and ST, and SW result. The symbols are retained after nondimensionalisation. 


SU 

mon2/?3 


+ Oj{^^S^ + + ^v^S-w' + 

•f — EAk\f EAc ak\0\v** co9$\ -f w^sin^i) 

EB\ - EBiS^jv** cos $i 4- w^^3in5i)j5^^ 

EC\^'* + EC2(v)^ cos $i — sin 


-F«a cos $i + ^Elg cos^ + Ely sin^ 

+ ^{Elg — EIy)w^*9m 2$i - £?Aca^v''cos^ 5l -f -w^^sin 25i^ 

+ (EAcak% - EB^y^' co9$i - F(72^"sin5ij5v'' 

+ sin $i + |f/s sin^ ^5 + + Ely cos^ 

+ ^(J5/z — £?/j/)v'^sin25i - EAe^^w" ain“ Si + ^v^sin25|^ 

+ (EAeak^ - ES^y /sin Si - EC 2 / cos Si 
4 sin 5i — w"cosSiJ 

4- (EIz - F/y)| "■ sin 2Si 4- v"ui'' cos 2Si | dac 


(2.17) 
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ST 




m||v + cos + 2y?pW + 2 jf (ijv^ •¥ ww') dx 
+ 2e|(t)^cQs^l 4* wSinifi) V + cp^8m#l]^o 


•« 

^p(* + 2v) 4* tw 4* cos 


Sw 


+ l{^m2 - + «y*(w^08#i - v'sin^i) 
Cposintfl 4* C(T/9pZ cos ^1 — Cp^i/sinli — tycos#!^!#^ 

- Cp^afcos^i 4- 2 iJco8^I^#w^ 

- Cp^zsin#! 4- 2vsm#i^5io'|</x 


(2J8) 


and 


SW 


i> vSv 4- Lu/Sw 4- M^S^) dz 


( 2 . 10 ) 


mon2i?3 JO 

The singly underlined terms are nonlinear and the doubly underlined terms are 
third order (c^) linear torsion terms. 


i 

t 

{ 

t 

j 

i 

i 

I 

i 
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3. AERODYNAMIC LOADS 

/ 

The external loads Lvi and of Eq. ( 2 , 10 ) are of aerodynamic origin, 
The expressions for these aerodynamic loads in hover, distributed along the length 
of the blade, are obtained based on a quasisteady strip-theory approximation. The 
flow is assumed to be incompressible and inviscid, The induced inflow is assumed 
to be uniform and steady. Virtual inertia effects are included. The variation in free 
stream velocity due to flap, lead-lag, and torsion motions are taken into account. 
Blade angle-of-attack stall is neglected. 

3.1 Circulatory Forces 

The forces L„ and are in the y and z directions respectively (i.e., in the 
undeformed coordinate frame). In the present analysis the strip theory is based on 
sections in the deformed frame. As a result, forces obtained in the deformed frame 
have to be transformed to the undeformed frame. The magnitude and direction 
of the relative wind over the blade section is obtained by vectorially adding the 
negative of the blade velocity to the induced inflow. 

+ (3.1) 

where V is the relative wind, % the blade velocity, and tJ,- the inflow velicity. Hodges 
and Dowell (1074) present the expression for P 5 in terms of the unit vectors in the 
undeformed frame. 

% ■■ (xi - Dj/i cosj9p)7x -f- [jii + (xi cos/9p - zi sinj^pjjtp 

+ (ii -♦•nyisin/3p)ij (3.2) 
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with 
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a?i ■■ at + « - f sinti) — w\tj e\n$i + f cos#]) 

VI ■■ V + 1; cos — f sin $i 

zi " w + i|sln#i + f cosl| 


(3.3) 


The following approximations are made. 

(а) is obtained at a reference point (t^nO) on the chord. 

(б) sin^p psf ^p] COS^p 1, 

(c) Terms of ord''T and higher are neglected. 

Then Eq. (3.1) is written as 

+ + (34) 

where 


C/* *“ ~ti -f l)^;r cos fii + w'tfr sin $i + n(w + tjr cos #i ) — v0p 

f/y ■“ —d + ^i/rsin^l — n(at + u — vVcos^l - w^t^psin^l) + 0^p{w + i^rsin^l) 

Ug =» — lb + cos 01 — ^^p{v + fjr cos ) ■" ‘'i 

We are interested in expressing V in the deformed frame since the airfoil section 
i>' considered in that frame. So if V" is written as 
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-Up 


U^\ 

-Ui> 

• - m 

uA 



U:) 


(3.7) 


The axial velocity Un is negligible in hover; the expressions for Ux and Up are 


Up * ” w^w^sintfi 

+ n[ucos^l + — n'/9p)co9^l -f («/^ + j9p)wsin^l] 

+ vcosli + (w + V|‘)sin^l (3.8) 


and 


Up “ —0*1^1 — + vWcostfi 

+ fl[-«sin^l - (vv'““ w^p)sin#| +(u/' + /Sp)(i7p + (i>co9#i)] 

« 

— 0 sin -f ( w + V,') cos #1 + 


(3.9) 


According to the quaaisteady approximation of the thin airfoil theory, the 
circulatory forces are based on the angle of attack at the three-quarter-chord point 
[t.g., see Fung (1969) pp. 193). Therefore U^ and Up are evaluated at the three 
quarter chord point; that is, tjr is the rj coordinate of the three-quarter-chord point. 

Figure 5 shows the aerodynamic environment in the deformed frame. The 
circulatory forces Zve, Iwo ^ directions, respectively, can 

be written as 
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Zwc « £# sin a — D cos a 
Xy^f. — « — Lcos a — D sin a 
— Muc 


(3.10) 


where L, and Mac lire, respectively, the lift, drag and pitching moment (about the 
aerodynamic center) per unit length, a is the angle of attack at the three-quarter- 
chord point, and is the aerodynamic center offset from the elastic axis. The 
expressions for L, D, and Mac nre 


L-.Ci^fV-c 

D - Cp\pVh (3,„) 

w.. - Cu„ipvh‘‘ 


with 


V-’^Uf + Up 
and 


tan a »= — 

Uj' 

where C/,, C/j, and are the section lift, drag, and pitching moment coefficients, 
and c is the section chord. Substituting Eqs. (3.11) in Eqs. (3.10) yields 
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-CdUtV) 

Iwc « -Ipc{ClUtV + CoUpV) 


♦ 

With the assumptions 


*= CO + cio 

Cd “ rfo + <^iOt + rfoa* 

C/p 

o #*sf sino ■■ ~— 
V/^Ut 


(3.12) 


(3.13) 


Eqs. (3.12) can be written as 


Zvc “■ + (co — di)UTUp + (ci — d2)Up 

Zwe " dij^UjUp + rf|l7^j (3.14) 

+ V^) - 


The assumed order of magnitude is unity for ci and d2\ t for cq, rfi, C'jv/ao 
<®/2 for and <“ for do- The circulatory forces L„c and Lu;c in the undeformed 
frame are obtained from Zt»c and Tjy,f. using the transformation matrix. 



fM 

- irl’’ 


.Zwe) 

^c, 
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(3.15) 


5 ." 

fl 

if 
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Thus, the circulatory aerodynamic force components to order are 


^ve 


-f(co - di)^— ^n*a* 8 in 2 tfi + Ha;t;|-cos 2 Ft^ 

+ (ci — rf 2 )|(Ox 8 inFi -*y,‘costfi)^ — 2 n“xi|rji?psin#ij| 

+ sin &i cos^ $i -f sin 2 tfi^ 

+ (ci + do)^Oxcos ^i(— Ox sin $i + V{ cos $i + ^ffrpp) — Oxvisin* $i -)- sin 2 Fi 

+rfl^O“x“ sin^^i — nxv;sin 2 ^l^| 

— un^x(do + 1/2) 2^1 sin $i 

- vfl^x^p[(ci — </2 )s*Q2(?icos $i + (ci + do)sin^l cos 2 Fi] 

+ «;n^a 5 ^p(do + 1/2) sin 2^1 sintfi 

+ «;'n^x»;,.(— Cl + do + 2d2) sin $i cos d\ 

v|cos(?i — 2 flxdo cos^ $i — flx(cQ — dj ) sin 29 \ +(c[ — d2)^20x sin^ ^1 — v,-sin 2 di^j 


+ sinFi|^ 20 xcocos^F| — (ci + doXOxsin 2^1 — Vfcos 2 Fi)+ 2 nxd[ sin^Fij j 
«;|costfl — flxdp sin 2 ^i + flx(co — di)cos 2 tfi +(ci — d2)^— Oxsin 2 tfi + 2 viCOs*^|^| 


+ sini?i Oxcosin 2 Fi +(ci + do)^Oxcos2^1 + v,-sin 2 Fi^— nxdisin 2 ^l|| 

+ ^Qxfjr{—ei + do + 2d2)sin0i cos^i 

— vv^Cl^x{do -H d2) sin $i sin 2 <?i 

— ott;^fl^x[— (ci — d2) sin 2 ^i cos $i -h (ci + dp) sin $i cos 20 i] 

— vro^a!^|(do + d2)sin^^lcostfi 




- 0 V 02*2 


(ci — d2)^sin2tfi cosidi ~ sin* tfi^ 


+ (ci + dp)(- »iB 2», CO. «1 +«■.»».)]) 


( 3 . 18 ) 
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i/»e^8intfi|-(/onVco9-^i +(co + n*v,’Cos2^i^ 

■F (ci - rf2)|(n*sintfi - v,-co9^i)“ - 2n’*«?/r/Ip8ln^ij| 

— costfi|co^n'^*"cos^(?l + nacv,‘8m2^i j 

-f*(ei + (/o)|naico8tfl(-nxsiii^l -f v»co8^i + Oijr^p) - Oxv; 8in“ #i + ~vjsin2tf| 
+rfl^n“*“ sin* ^1 — flafv,* sin 2^i 

- u2n**J(ci ~ rf2)sin*tf| +(ci + rfo)sin^l cos'^j 

— t>n“*i^p((ci - (/2 ) sin sin 2^1 4-(ci + (/o)co8tfi co8 2tfl) 

- w 2n"a;?p[(ci “ t/2)sin® + (ci + rfo)sin#l co8“ 0i 


— u}'^'‘xrir 


2(ci - rf2)8in" + (ci -h dQ)cos“ 9i 

t ■* . 

-F t)|8intfi|^ — 2nz(/ocos^ flx(cQ — <f|)9in2tfj +(ci —• rf2)^2fl*sin"^i v,‘sin2tfi^| 

- C08 ^l|2na:co COS' -(ci +t/oXOx8in 2^1 “V,- cos 2^l) + 2n*rfi sin'^ijj 

■f w|sin<^i^ — flxrfosjn2^i + nz(co — <^i)cos2tfl + (ci ~ rf2)^— Oz sin2»,+2uiCOj=«i)] 

— cos^ifnzcosmS^t + (ci +t/o)^na:cos2tfi + tt,‘sin2^l^— Oztii sin2tfi|| 

2(ci - </2)sin“ $i + (ci 4- rfo)cos“ 


— ^Clxijr 

- vv'n^z 


2(«i — <f2)sin®tf| +(ci + do) sin 2^1 cos^i 


— nw^n“z((ci — d2)sin^t sin 2^1 4- (ci 4-do)cosi^i cos2^l] 


-o'‘02^2 


(c| — d2)sin®^l 4-(ci 4- do) sin cos' 9i 


4- o^w^fl'z^[(ci -• d2)sin^l sin 2^1 4-(cj 4- do)cos^i cos2(?i] 


<2 2 •> 1 

— n^z*~ 

m 


(ci — d2)sin® ^1 4- (ei 4- do)sin $i cos 


!d> 


(3.17) 
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i/»c^n-ac- jcCA/ac “ **0) sin cos | - d| sin^ |j j | 

- Uxeji{ei f (/o)^t)sin2^i -wco92$^ 


( 3 . 18 ) 


The singly underlined terms are nonlinear and the doubly underlined terms are 
third order damping terms. 

3.2 Noncirculatory Forces 

The noncirculatory components are obtained from unsteady thin airfoil theory 
(Fung (1960), pp 209>210]. To order <*, they are 

f 1 2 

Lxunc 

- \^P<^“ 


-ib + rSac^ ( J* j^J 

Qe + ejV - Qc -I- 1 


( 3 . 19 ) 


3,3 Combined Forces 

The expressins for the complete aerodynamic forces in the lead-lag, flap, and 
torsion directions are given by 


Lv ^ Lv<; 

"F f^ne 
A/^ " "I” 


( 3 . 20 ) 


In the derivation of the above aerodynamic forces, nonlinear rate product terms 
such as v^, if, t/w are neglected while retaining other nonlinear terms up to order 
e“. This is necessary for a linearized perturbation analysis about the steady state. 
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Some damping terms or order which may bo important for the stability analysis, 
are retained. The uniform induced velocity is based on the momentum theory. 
The pitch of the blade is obtained by combining the momentum and blade element 
theories. 



•i 


where kh is an empircal factor, Cf the thrust coeflicient, a the solidity ratio of the 
rotor, and ^75/; the blade pitch at three*quartcr span. 


6 
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4. FINITE ELEMENT DISCRETIZATION 


4,1 Equations of Motion in Terms of Nodal Degrees of Freedom 

The finite element analysis is used to discretize the spatial dependence of the 
equations of motion. The blade is divided inJo a number of beam elements. Every 
element consists of two nodes, denoted by node-1 and node-2 (Fig. 6), with five 
degrees of freedom namely v, v', w, w', and ^ at each node. The rc][^on for choosing 
^ over ^ as the torsional degree of freedom is as follows. The aerodynamic loads 
depend explicitly on If is chosen os the torsional degree of freedom, then 
the presence of the integral term /q* v^w*dz [Eq. (2.2)] makes the global matrices 
nonbanded because of the coupling of the degrees of freedom of one clement with 
that of the other elements. The choice of ^ eliminate this integral and thus 
preserves the banded structure of the global matrices. 

Hamilton’s principle, Eq. (2.12), is discretized as 


(4.1) 

1-1 

with 

Aim6Ui--6Ti-SWi (4.2) 

where SUi, 6Ti, and are respectively the strain energy, kinetic energy, and the 
virtual work contributions of the ith element and n is the total number of elements. 
The distribution of the deflections o, w, ^ over an element is represented in terms 
of the nodal displacements using shape functions. For the ith element (Fig. 6), 


- ee 
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PI M 


(4.3) 


where the shape function matrix p] is 



Hi Ho «3 0 0 0 0 0 O’ 

0 0 0 0 Hi Ho H 3 H4 0 0 

00000000 ^^i 


and the vector of element degrees of freedom {g»} is defined as 


(4.4) 


{9») “ 1*^1 *^l ^ <^l ^l ”^2 h 

The nodal degrees of freedom at node-1 of the element are vi, wi, w\, h and 
those at node-2 axe vo, W 2 , wL, ^o. The shape functions in Eq, (4.4) are the 

Hermite polynomials defined as 



(4.6) 


H^(®i) 
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I 


where is the length of the ith element and *,• is the local axial coordinate for 
the ith element, measured from the left end of the element (Fig. 6). Similar to 
Eq. (4.3), the distribution of the virtual displacements Sw, and over the ith 
element is assumed to be 






PI (»«.•> 


(4.7) 


The substitution of Eq. (4.3) and Eq. (4,7) in Eq. (4,2) results in 




(4.8) 


which is nonlinear in qj. Thb function can be conveniently written as 




(4.9) 


where [Mdq,)j, (C,-(q,)l, and [K',(7,)), represent the element inertia, damping and 
stiffness matrices, respectively and {(2,} is the element load vector for the ith 
element. The global matrices are obtained by the assembly of the element matrices. 
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The global degree of freedom vector is denoted by {q) and the global virtual 
displacement vector by {^ 7 }. The global load vector (Q) is formed by combining 
the element load vectors {Q,}. The assembly of element inertia matrices (Mdijj)] 
results in global inertia matrix (A/(< 7 )j as shown below. 


[M(q)l 



The matrix [M] is a square matrix of order 5(n + 1 ); it is also banded with the 
semibandwidth N, equal to 10. The global damping and stifness matrices (C'( 9 )] and 
[/C(g)] are obtained in the same manner. The matrix [/C] also is banded, whereas 
the matrix [C] is nonbanded because of the presence of the double integral term 

lo lo ^ dsc dxi 

in the expression for 6Ti [Eq. (2.18)]. Matrix (M] is symmetric; the presence of 
terms in [C| and [K] due to the nonconservative aerodynamic forces makes them 
asymmetric. The process of assembling element matrices to obtain global matrices 
can be mathematically indicated by substituting Eq. (4,9) in Eq. (4,1), 
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A - {{«)'’|Af(,)l(«} + {f«)'’|cr.(,)lW) 

+ {<,)''| K (,)]{,) -0 


( 4 . 10 ) 


Since the virtual dispalcements Sq are arbitrary, Eq. (4.10) leads to the equations 
of motion 


[M{q)]{q) + {C{q)]{4) + [K{q)]{q) - {Q) (4.11) 

These equations are nonlinear in q. The bandedness of (A/) and [/C] helps reduce 
the storage space needed in computer programs since the zero elements off the band 
need not be stored. 

4.2 Boundary Conditions 

The formulation of the problem is based on energy pricniples; and hence, the 
force boundary conditions, which are imbedded m the formulation, do not have to 
be considered separately. The nature of the displacement boundary conditions at 
the root determines the configuration of the rotor. For example, v, w, «/, and ^ 
are zero at the root for a hingeless rotor blade; and for an articulated blade, v, w, 
and 4 are zero at the hinge. 
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5. SOLUTION PROCEDURE 


5.1,.Ttim-.SQliitiQa 

The steady trim condition of the rotor in hover is determined first since the 
flutter stability equations are linearized about the trim state. The nonlinear steady- 
state equations are obtained by dropping all time-dependent terms from Eq. (4.11). 


(li^o(9o)I{go} “ {<3o} (5.1) 

The matrix [Kq\ is asymmetric, banded, and a function of the steady displacement 
go- The contribution to the matrices [iirol {^o} is from both structural and 
aerodynamic forces. The expressions for the terms of the steady-state element 
matrices are given in Appendix C. These matrices are evaluated numerically us- 
ing Gauss quadrature formulas. The numerical solution of the nonlinear steady- 
state equations is evaluated iteratively using Brown’s algorithm [Brown and Dennis 
(1972)). This algorithm is a modified form of the Newton-Raphson method and the 
solution is estimated such that the sum of the squares of the errors in the equations 
is a minimum. The linear solution of Eqs. (5.1) is used as the initial estimate in 
the solution procedure. The nonlinear stiffness matrix [/Co(go)l is updated at every 
step of the iteration using the estimate from the previous step. The bandedness of 
[itol is useful in two ways. The storage space needed in the computer program is 
less since only the elements within the band are stored. The number of computer 
operations is reduced considerably since the matrix operations corresponding to the 
zero entries of [iCob outside the band, are skipped, 

For a given level of thrust, the collective pitch required at three-quarter span 
(^.75i?) is calculated from Eqs. (3.18). The steady elastic twist modifies the pitch 


- 31 - 


distribution of the blade. The new pitch at three-quarter span is *^.75/? + ^,75/? and 
corresponding to this pitch vallue, a new thrust level is calculated. 

5.2 Flutter Stability 

The flutter equations of motion are linearized about the steady-state solution. 
The displacement q is written as the sum of a steady component {qo) and an unsteady 
perturbation (9). 


q>^qo+q (5.2) 

Substituting Eq. (5.2) into Eqs. (4.11), subtracting Eqs. (5.1), and keeping only 
linear perturbation terms, the flutter equations are obtained as 


A/(<7o) I?! ■+• C{qo) + I<{%) 


{9} “= {0} 


(5.3) 


The inertia, damping, and stiffness matrices are functions of the steady deflection 
qo. The expressions for the terms of the perturbation element matrices are given in 
Appendix D. 

The flutter stability is determined by studying the eigenvalues of Eqs. (5.3). 
The normal mode method (which considerably reduces the number of equations) is 
used to calculate the flutter eivenvalues; and hence the natural modes are needed. 
The coupled rotating vibration characteristics about the equilibrium position are 
evaluated by dropping the damping matrix and removing all aerodynamic terms 
from Eqs. (5.3). The resulting equations are 


=» {0} 


(5.4) 
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The inertia and stiffness matrices [A/,] and |/sT,) arc symmetric and the eigenvalues 
(frequencies) of Eqs. (5.4) are real. 

The flutter equations of motion, Eqs. (5.3), are transformed to the modal space 
by writing 


{q} - (4>]{p} (5.5) 

where (<!>] is the matrix of the first N eigenvectors (colulmnwise) and (p) is the 
vector of N generalized coordinates in the modal space. Substituting Eqs. (5.5) 
into Eqs. (5.3) and premiiltiplying by [4»]^ yields the modal-space equations 


(A/*Kii) + Ic'Hp} + [ic*|(p) - (0) (S.0) 


The matrices (Af*l, [C7*], and (iC*) are asymmetric and of order N. These 
equations can be written as a first order system of 2N equations. 


[A]{r) - iBl(r) - (0) 


where 




w m 

. (Oj (/I . 


■(0! [K*\ 

.(/) 10 ). 


(5.7) 
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By writing r as 



(5.8) 


Eqs. (5.7) becomes 


!Al(r)-X(BI(f) (5.9) 

These equations are solved as an algebraic eigenvalue problem. The eigenvalues X 
are complex. 



0. APPLICATION TO SINGLE-LOAD-PATH BLADES 


To show the fesibility of the finite element formulation, the method is first 
applied to single-load-path blades; a hingeless blade and an articulated blade with 
a hinge offset of 5%are considered. Convergence studies with respect to increasing 
number of finite efiements and natural modes are presented. For the hiugeless blade, 
comparisons are made with previous theoretical results obtained by modal methods 
and with some experimental results. 

The numerical results obtained are for blades with uniform spanwise poperties. 
The chordwise offsets of the center of mass, aerodynamic center, and tension center 
from the elastic axis are considered to be zero. The section constatnts EBi, EB2, 
and the warping constants ECi, EC 2 , are taken to be zero. A precone (^p) of 0.05 
rad ( 2.0 deg). Lock number ( 7 ) of 5, and solidity ratio (<r) of 0.1 are used. 

6.1 Hingeless Rotor Blades 

The uniform blade properties selected for the stability analysis of a hingeless 
rotor blade are given in Table 1. The stiffnesses Ely, EIz, GJ, and the inertial 
parameters chosen such that the rotating frequencies corresponding 

to given values. The rotating flap frequency of the blade is taken to be 1.150. Two 
different lead-lag frequencies are considered; a frequency of 0.70 represents a soft- 
inplane blade whereas 1.50 represents a stiff-inplane blade. Similarly a torsional 
rotating frequency of 2.50 represents a torsionally-soft blade and 50 a torsionally- 
stiff blade. 

The convergence of the steady state deflections using different numbers of finite 
elements is presented first. Table 2 shows the steady tip deflections 
(nondimensionalized with respect to the rotor length) and (in rad) as the total 
number of finite elements is varied from 2 to 8 . The results are for a thrust level, 
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(7r/<T, of 0.1. The results show, for the case considered, that six elements are 
sufficient for good convergence (three to four digit accuracy). The banded nature of 
the stiffness matrix reduces both the storage space and computation time required. 
For example, if eight elements are used, the full storage mode would require 2025 
storage locations whereas the band storage mode needs only 855 locations for the 
stiffness matrix. Furthermore, the number of arithmetic operations involving the 
stiffness matrix is reduced by about 60%. The use of the linear solution as the 
starting vector makes the nonlinear iterative solution converge rapidly; only three 
or four iterations are needed for a converged solution. As a result, the computation 
time required for the trim solution is small. 

Figure 7 shows the steady tip deflections of a stiff-inplane hingeless blade as the 
thrust (and hence the pitch is varied. The torsion frequency corresponds 
to a torsionally-stiff blade for these results. The magnitue of the flap and lead-lag 
tip deflections increase with thrust. The magnitude of the torsional tip deflection 
increases first and at higher thrust levels decreases. The torsional deflection is 
primarily due to nonlinear flap-torsion and lag-torsion couplings. The steady-state 
aerodynamic torsional moment is zero since and are assumed to be zero 
[see Eq. (3.18)). At low values of Cy/ff, both the flap and lead-lag deflections are 
negative. At high values, the flap deflection is positive while the lead-lag deflection 
is negative. The results of Hodges and Ormiston (1976) are also shown in this figure. 
Hodges and Ormiston use a modal method with five nonrotating beam modes for 
each one of the deflections wo, wo, and ^o- The agreement between the two results 
is excellent except at high trust levels where a slight deviation appears. This may 
be due to the assumption by Hodges and Ormiston that sin0i^$ and cos(? 1 (not 
assumed here) for the claculation of aerodynamic loads, which may not be very 
accurate at high values of $. The steady tip deflections with respect to thrust for 
a soft-inplane blade are shown in Fig. 8. The flap, lead-lag, and torsion curves 
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are similar to that of the previous case. The flap deflections are almost identical 
to that of the stiff-in plane blade. The lead-lad deflections are higher in magnitude, 
as expected, since the lead-lag stiffness is lower. The torsional deflections also are 
higher in magnitude than for the stiff-inplane blade, which can be explained os due 
to the coupling between lead-lag and torsion modes. Corresponding results from 
Hodges and Ormiston (1976) are included in Fig. 8. Comparisons show that, as 
before, the agreement between the two results is very good at low thrust levels and 
the deviation increases with the thrust, 

The rotating coupled natural frequencies of the stiff-inplane blade about its 
steady deflected position are investigated for convergence as the number of finite 
elements is varied. Table 3 presents the fundamental lead-lag, flap, and torsion 
frequencies with increasing number of elements. It is seen that six elements are 
sufficient for four digit accuracy. The bandedness of the inertia and stiffness 
matrices is used effectively to mimimize computation time. 

The solution of the flutter equations results in complex eigenvalues. The normal 
mode method is examined for convergence as the number of modes is varied, keeping 
the number of elements fixed at six. The real parts of the first three eigenvalues for 
different numbers of modes is given in Table 4, This table shows that five modes 
result in a well-converged solution (five-digit accuracy). Figure 9 shows the root 
locus plot of the fundamental lead-lag, flap, and torsion modes as is varied 
from 0 to 0.3. Even though blade angle-of-attack stall has been neglected, such 
high thrust levels are included for comparision with other theoretical results. Six 
finite elements and five normal modes are used in obtaining these results. The flap 
and torsion modes are stable over the entire range of C-j'l<r considered, whereas the 
lead-lag mode is unstable for Cy/tr between 0.01 and 0.05 and above 0.17. The 
complete solution using six elements and five normal modes at a particular thrust 
level requires about 7 s of computation time on a CDC 7600 machine. Figure 
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10 shows the flutter stability boundaries as the lead-lag frequency is varied. The 
results of Hodges and ormiston (1070) are also shown for comparison in this figure. 
In this reference, the stability results are calculated using the normal mode method 
with five coupled rotating modes. There is general agreement between the two 
results. The differences that appear between the two results may be due to change 
in aerodynamic forces by including all second-order terms (order e-) in the present 
analysis. 

Figure 11 shows the root locus of the lead-lag mode of a stiff-inplane hingeless 
blade when the blade pitch is varied from -10 deg to 12 deg. The blade precone 
is zero for this case. Experimental results from Sharpe (1983) are also shown in 
this figure. In the theoretical analysis structural damping is included such that the 
results match that of the experiment at zero pitch. It is seen that the agreement 
between the present analysis and experiment is good at low pitch settings and not so 
good at higher pitch settings. Other theoretical analyses show similar discrepencies. 

6.2 Articulated Blade 

The same blade properties as given in Table 1 are used for the stability analysis 
of an articulated blade with a hinge oflset of 6%. The convergence characteristics 
for the nonlinear steady solution, natural frequencies, and flutter eigenvalues are 
presented in Tables 5-7, respectively. These results confirm the conclusions reached 
for the hingeless blade that about eight elements and five normal modes are adequate 
for four-digit accuracy. As before, only a few iterations result in a converged 
nonlinear trim solution. The equilibrium tip deflections are plotted with respect 
to thrust in Fig. 12. The lead-lag tip deflections are much higher than that for 
the corresponding hingeivess rotor blade because the lead-lag stiffness is lower when 
the blade is hinged. The root locus of the flutter eigenvalues obtained using six 
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dements and five normal modes is shown in Fig. 13. All three modes are stable 
over the range of Cj/ff shown. 



- 89 - 


7. MODIFICATIONS FOR ANALYSIS OF MULTIPLE-LOAD-PATH BLADES 

Having shown that the finite element formulation can be successfully applied for 
the determination of the aeroelastic stability of single-load-path blades, the analysis 
is extened to include the more complex multiple-load-path blades. As indicated in 
the Introduction, the flexure of a bearingless blade contains one or two flexbeams 
and a torque tube. The analytical model of a bearingless blade is shown in Fig. 14. 
Each of the flexbeams and the torque tube are modelled as individual beams. For 
this purpose several modification have to be made to the formulation presented in 
the previous chapters. 

7. 1 Inclusion of Axial Degree of Freedom 

In the analysis for single-load-path blades, the axial deflection u has been 
eliminated in terms of the other deflections and the expression for the centrifugal 
force. The force equilibrium equation in the axial direction cannot be solved a priori 
for multipie-load-path blades; and hence the u deflection has to be included as a 
degree of freedom. 

7.2 Refined Finite Element 

The beam element considered for discretizing multiple-load-path blades is shown 
in Fig. 15. This element is more refined than the simple beam element used for 
single-load-path blades. The element includes the axial degrees of freedom. Each 
element consists of two end nodes and three internal nodes with a total of fifteen 
degree of freedom. Each of the end nodes (1 and 2) have six degrees of freedom, 
namely u, v, i/, w, «/, and There are two internal nodes for u and one inter- 
nal node for Two internal nodes are needed for the axial degree of freedom 
to accurately represent the centrifugal force distribution over the element. With 
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two internal nodes, the distribution of u over the element is cubic. Hence, the axial 
force distributiion is quadratic and of the same order as that of the centrifugal force. 
The internal node for the torsional deflection assures that the torque approxima- 
tion over the element is of the same order as the flap and lead-lag bending moment 
approximations. 

7.3 Displacement Compatibility Conditions at the Clevis 

The element matrices are assembled, by considering displacement compatibility 
conditions at interelement boundaries, to form global matrices. Except at the clevis, 
the displacement compatibility conditions generally imply the continuity of u, v, 
w, w*, and ^ across element boundaries. The root element of the outboard blade is 
connected to the outermost elements of the flexbeams and the torque tube through 
the clevis which is assumed to be rigid (Fig. 14). The inboard beams are numbered 
1,2,3, etc. The displacement compatibility conditions at the clevb can be written 
as 

u; =» u — rjiv* 

Vi ■* V 

f7.11 

Wi^mw + iJi{$ae + ^) 

W*i ■■ 

and 

^i “ ^ac + {forfUzhtama) 

^ {fortorquetubt) (7.2) 
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where Vf, and are the deflections at the clevis for the ith inboard beam; and 
u, 4 , w, and <f> are the deflections at the clevis for the outboard beam. The 17, • is the 
elastic axis offset of the ith inboard beam from that of the outboard beam (positive 
forward) and Sac the difference between the collective pitch and the pretwist at the 
clevis (Fig. 16), Equation (7.2) for the flexbeam indicates that a rigi<l body pitch 
.rotation of the outboard blade is transferred to the flexbeam as an elastic twist. 

7.4 Modifications in Solution Procedure 

In addition to the differences noted above, there are two major differences in 
tL <;ion procedure for nonlinear trim solution. Some of the terms of the matrix 
(/iTol depend on the centrifugal force F, The centrifugal force distribution is known 
a priori only over the onboard blade; it is unknown over the flexure because of 
the multiple load paths involved. In obtaining the linear estimate of the steady 
equations, it is assumed that the centrifugal force in each of the inboard beam is 
in the ratio of its tensile stiffness EA. Then, at successive iterations the centrifugal 
force distrubution is updated using Eq. (2.6). 

The lead-lag and flap stiffnesses depend on the orientation (pitch) of the cross 
section. The pitch distribution is known a priori for the outboard blade. For the 
flexbeams, the resultant pitch is the sum of the pretwist 9pc and the elastic twist due 
to 9ae at the clevis (Fig. 16). For the first iteration of the nonlinear trim solution, 
the pitch distribution along the flexbeams is assumed to be linear. This distribution 
is updated at successive iterations based on the torsional deflection. 


8. APPLICATION TO MULTIPLE-LOAD-PATH BLADES 

The formulation modified for multiple-load-path blades is applied to several 
analytical models of bearingless rotor blades in hover. The outboard blade, the 
flcxbeams, and the torque tube are considered to be uniform, The chordwise offsets 
of the tension center, mass Center, and aerodynamic center from the elastic axis are 
assumed to be zero. The section constants EBi, EB 2 , and the warping constants 
ECi, ECo are taken to be zero. The aerodynamic forces over the flexure is neglected, 
A prccone (/3p) of 0.05 rad (2,9 deg). Lock number (7) of 5, and solidity ratio (<r) of 
O.I are used. The offset of the elastic axis of the fiexbeams (Fig. 14) from that of the 
outboard blade (i» »■ -tfz) is taken to be 0.4c, where c is the chord of the outboarad 
blade; the torque tube offset (>72) is zero. The length of the flexure is 0.25/?. The 
blade structural poperties are given in Table 8; the aerodynamic properties are the 
same as given in Table 1. Stiff- and soft-inplane and torsionally-stiff and soft blades 
are considered. 

gil-T w ip be am M odel wit h Z m. Jnbi?.ar.d--Pi.tfih 

To understand the behavior of multiple-load-path blades, a simple analytical 
model is considered first. This model has an inboard flexure cosisting of two identical 
flexbeams; the torque tube is not included (Fig. 17). The pitch of the flexure is 
fixed at zero while the pitch of the outboard blade is varied. This is equivalent to 
a physical model with a pitch bearing at the clevis. First, a convergence study is 
made to determine the number of finite elements needed to obtain the nonlinear 
trim solution to a desired accuracy. Table 9 shows the stead> state tip deflections 
vojjp, (nondimensionalized with respect to the rotor length), and 
rad) as the number of both the flexbeam and the outboard elements is varied. The 
blade considered is stiff inplane and soft in torsion. A thrust level, of o.l is 
used. The results show that two elements each for the flexbeams and two elements 
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for the outboard blade are sufRcient for five digit accuracy. The nonlinear iterative 
solution converges rapidly since the linear solution is used os the initial estimate; 
only four or five iterations are needed for a converged solution. 

The rotating coupled natural frequencies of the blade about its steady defleetd 
position are inve^itigated for convergence as the number of finite elements is varied. 
Table 10 presents the fundamental lead-lag, flap, and torsion frequencies with 
increasing number of elements. Again, it is seen that a total of six elements, two 
for the outboard blade and two each for the flexbeams, are sufficient for Ave digit 
accuracy. 

The solution of the flutter equations results in complex eigenvalues, The normal 
mode method is examined for convergence as the number of modes is varied, keeping 
the total number of elements fixed at six. The real parts of the first three eigenvalues 
for different numbers of modes is presented in Table 11. This table shows that five 
modes result in a converged solution (five digit accuracy). 

Figure 18 shows the root locus plot of the fundamental lead-lag mode as Cj'/tr 
is varied from 0 to 0.2. It is seen that the lead-lag mode is unstable for CTj<r higher 
than about 0.05. In this figure, the results from an equivalent-beam model are also 
shown. In general, the properties of the equivalent-beam model ore obtained by 
matching the fundamental frequencies of the equivalent-beam model with that of 
the multibeam blade. For flexbeams with uniform spanwise poperties, analytical 
expressions for equivalent properties can be obtained. These are given in Appendix 
D, The two results shown in Fig. 18 diflTer from each other, particularly at low 
and high levels of thrust. The results at low thrust levels are of significance to tail 
rotors. A careful study has been made to determine the cause of the differences 
between the two results. The equivalent-beam model does not accuravcly simulate 
the cross-coupling stiffness terms which are nonlinear in nature and depend on the 
steady deflections. Results have been computed by suppressing, in the perturbation 



analysis, nonlinear cross-coupling terms of nonlinear origin from the structural 
stiffness matrix of the flexure, for the twin- beam as well as the equivalent beam 
model. These are shown in Fig. 10. Both sets of results are identical. 

For the same model but considering a soft-inplane, torsionally-stiff blade, the 
root locus of the lead-lag mode is plotted in Fig. 20, The behaviur of this blade 
is very different from that of the stiff-inplane blade. The coressponding equivalent- 
beam model results are also included in this figure. The discrepency between the 
results of the two models are more pronounced, particularly at high thrust levels, 
than that of the stiff-inplane blade. Again the removal, in the perturbation analysis, 
of the cross-coupling terms of nonlinear origin from the structural stiffness matrix 
of the flexure makes the results identical (Fig. 21). Similar conclusions are reached 
from the root locus plots of the lead-lag mode of a torsionally-stiff blade (Figs. 22 
and 23). 

8.2 Twinbeam Model with Same Pitch For Inboard and Outboard Segments 

This blade model also has a flexure consisting of two identical flexbeams; and 
there is no torque tube. Both the inboard and the onboard segments are maintained 
at the same uniform pitch distribution. This represents a physical model with pitch 
bearings at the hub end of the flexbeams. The root loci of the lead-lag mode with 
respect to increasing level of thrust are shown in Fig. 24 for both the twinbeam and 
the equivalent-beam models. The results are for a stiff-inplane and torsionally-soft 
blade. The two results are in major disagreement with each other. This disparity 
is explained as follows. The fundamental frequencies of the equivalent-beam model 
are matched with that of the twinbeam model at a chosen pitch setting (zero in the 
present example). As the pitch is changed, the effective flap and lead-lag stiffnesses 
vary because of flap-lag coupling. This modifies the frequencies of the blade; and 
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the frequency match no longer exists. This phenomenon is clearlyl seen from Table 
12 where the fundamental frequencies are presented for various tbjrust levels. As 
the thrust (and hence the pitch) increases, the difference between the frequencies 
of the two model widens. Thus, the equivalent-beam simulation is not accurate 
with changing pitch. The results of Fig. 25 are for a soft-inplane blade. Again the 
discrepency between the twin-beam and the equivalent-beam models is big. Similar 
effects are observed from the results of a torsionally-stiff blade (Fig. 26). 

8.3 Bearingless Model with Torque Tube 

This blade has a root flexure consisting of two identical flexbeams and a 
torque tube in the middle. Thb model is a more appropriate representation of 
the bearingless blade. Pitch control of the blade is achieved via the torque tube 
by rotating the clevis. Due to the piltch application at the clevis, there is an 
elastic pitch distribution over the flexbeam and this is determined iteratively. As 
a result, the number of iterations required tor convergence of the nonlinear trim 
solution is increased. Now it takes about ten iterations for five digit accuracy, 
whereas, four or five itereations have been adequate for the two previous models. 
The convergence characteristics are similar to the one presented for the previous 
models in tables 9-11. The results have indicated that a total of eight elements - 
two elements for the outboard blade, two for the torque tube, and two each for 
the flexbeams - are adequate for five digit accuracy. The root loci of the flutter 
eigenvalues (fundamental lead-lag, flap, and torsion modes) are shown in Fig. 27 
as Cj'ItT is varied from to 0.15. The flap and torsion modes are stable while the 
lead-lag mode becomes unstable at high levels of thrust (Cj/^ greater than about 
0.13). This is a weak instability and the introduction of a small amount of structural 
damping would stabilize the mode. 
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0. GENERAL CONCLUSIONS 


The finite element method has been successfully applied to determine the 
nonlinear trim deflections and the flutter stability of sing! v and multiple-load>path 
blades in flap bending, lead-lag bending and torsion. The formulation is based on 
Hamilton's principle, The spatial dependence of the equations of motion is made 
discrete by dividing the blade into a number of beam elements. Nonlinear trim 
deflections are evaluated iteratively, solving the complete set of of global equations, 
without making a modal transformation. The solution procedure is made efficient 
using the bandedness of the stiffness matrix and the linear solution as the starting 
vector. The coupled rotating natural modes are calculated about the trim condition 
of the blade. The normal mode method based or. these modes is used to solve the 
linearized flutter equations as an eigenvalue problem. 

The formulation is first applied to uniform single-load-path blades, namely 
hingeless and articulated. Six finite elements are sufficient to determine the trim 
deflections, and the free vibration characteristics. Five coupled rotating modes 
are adequate for the normal mode solution of the flutter equations. The steady 
tip deflections agree very well with that of a previous research which used modal 
method based on a total of fifteen nonrotating uncoupled modes; the flutter stability 
boundaries also show axcellent agreement with that previous research based on a 
modal approach. The root-locus plots of flutter eigenvalues as a function of the 
pitch are generally in good agreement with experimental results. 

Suitable modifications have been made in the formulation so as to consider 
multiple-load-path blades. Numerical results are presented for several baringless 
models. Results of the redundantbeam analysis are compared with that of the 
equivalent-beam modelling; and the following important conclusions are reached. 
The equivalent beam does not simulate the nonlinear structural coupling terms ap- 
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propriately. Also, the properties of the equivalent beam are estimated by matching 
frequencies with the bearingless blade at a particular pitch setting. At other pitch 
settings, the frequencies change and the match no longer holds. As a result, the 
stability results obtained using the equivalent-beam model are inaccurate. 

Future extensions may include forward flight, ground and air resonance, and 
consideration of composite blades. 


a 
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APPENDIX A, NONDIMENSIONAL QUANTITIES 


The nondimensional quantities in SU, ST, and 6W [Eqs. (2.17), 
(2.18), and (2.19)) are defined here. These symbols are retained 
after nondimensionalisation. 


Old quantity 

New quantity 

Assumed order 
of magnitude 

xfR 

X 

1 

m 

() 

1 

vjR 

V 

c 

w/R 

w 

c 



c 

9 

9 

1 

F/mon^R" 

F 

1 

Elylmon-R* 

EI,t 

1 

Elflmod'R* 

El, 

1 

GJlmon^'R* 

GJ 

1 

EAfmon^R^ 

EA 

1/e* 

EBilmod^R^ 

EBi 

ft 

c* 

EBsImoCl-R^ 

EBz 

e 

ECilmon^R^ 

ECx 

e* 

ECifmoil^R^ 

ECz 

e 

K'aIR 

kA 

e 

^m\fR 


e 

f^mtlR 


e 

ejR 


,8/2 

c,IR 


,8/2 

P* 


e 
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APPENDIX B. STEADY-STATE ELEMENT MATRICES 


The expressions for steady-state element stiffness matrix and load vector are 
derived from the steady part of Eq. (4.2). 


Ao.mSUo.-’STo.^SWo. 


(BA) 


The strain energy, kinetic evergy, and the virtual work contributions of the ith 
element in steady state are obtained from Eqs. (2.17), (2.18), and (2.10) by dropping 
all time dependent terms and writing 
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ORIGINAL RAGE 18 
OF POOR QUALITY 
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where the steady-state aerodynamic forces are 
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+ voivox®((ci — t/ 2 )sin^sin 2 tf + (ci + rfo)costf cos2tf| 

1® 2 1 
- tvS'x®- 

(ci — rf 2 )*in® 9 + (ci + (fo)sin tfcos® 9 

) 


(B.7) 
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’^*0 ~ 5*°{‘ 




(ci + </o)»in^cos# 


- d\ sin" 


{B.9) 


The Ao,- of Eq. (B.l) can also be written as 


Aof “ {%}rao(9.)o + (^‘®) 

By comparing Eqs. (B.l) and (B.O) the element stiffness matrix and the load vector 
can be written. The contribution from 8Uq. and 6Tq^ to the element stiffenss matrix 

and the load vector is given as follows (the contribution from ^VVb,. can be easily 

inferred from Lt»o» 



■ 






K%ov 

[iCtvwi 


• 


[*:*„) 



and 


(B.IO) 


where 


({Qv}] 

(Qilo ■** { 

U«*)J 


(B.ll) 




ffolB') {H')^ + {Eh cos® < + El, siB® e - EAel cos »® j{H") {«*)’' 
m{B}(W>®'j iix( 


\Kvw\ 



EAtfj sin 2${H"] {H^fdxi 


- 55 ~ 



ORIGINAl PAGE IS 
OF POOR QUALITY 

+ Fo«a8in<{^^}{Hfr 

-EC's sin 

+ mej,flin^(JY}{//^}^ 

— m«^af8ln^(/ir^}{//^}^ 

- (Eh - /?/j,)sin2^{H''}{H"}^{vS){i/^}^ 
+ (E/, - EIy)tos2d(H^')[H>>f(v}l)[H^f 

+ cfy{//»){fl')''{«’!>{«i}n<'». 


|iCu;v! * l'^««/l 


[iCtt^w] - ’ ^Fo{H'){H^f +(^Eh sin- 9 + Ely cos^ $ - EAel sin^ rfx.- 


[iCu,^] - ]^EAeak\ - £;5o)<>' 9in 6{H''] 
-FoeaCos$[H^^](R^f 
+ ECo C09 

+ mefXC09$(lF}(H^}^ 
■^■(Eh-EIy)co92${H'*}{H^f{vi}{H^f 
+ (Eh - EIy)sin20{H^}{H'*f(w*Q){H^f 

+ o/(h'){h"}’'{«5}|hv}’^ dxi 
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-f m«j sin ${H^}{H)^ 

— megX 9 \n 9 {H^}{JHl*}^ 

- i(B/, - BIy)«nU{H^){vif {»»]{«'< f 

+ (Eh - dzi 

lK-^1 - jf'* \{EAo,k^ - BBiy sin 

+ E<72Cos»|h}|{H*)^ 

+ tnefXco 59 {H^){H')'^ 

- i(B/. -B/,)sin2»{ff^}{«5}^{H*){//')’’ 

\K^ 9 \ - fl' f( 0 '^ + - EAk*y + 

+ EGi|«{||h}|’^ 

+ w(*m2-*m,)«<»2<(W^){H^)’^ d^i 

tU 

{Qv} ■» / [Fo«o cos^{iif^} + mtg costf{H} — megXcos9{H*)] dx; 

JO 

fU 

{Q'w} ■" / [FQea»^f^9{H^] + m^pX{H] - fntgX9\n9[H^]] dx{ 

JO 

{Q<f,] “ “ - rntg^pX cos 9 {H 4 ,}^ dx{ 
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The column vectors (H), {H^}, and (ej) are defined as 

{»)*■- Iffl Hi H) Hi\ 

{e})*' — l«o, ve, “Oj vJjJ 
A similar definition holds good for {«;$} . 
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APPENDIX C. PERTURBATION ELEMENT MATRICES 


The strain and kinetic ener^ contributions to the perturbation inertia, damp- 
ing, and stiffness element matrices are given here, The contribution of the aerodynamic 
forces is easily inferred from Eqs.(3.16)-(3,18). Since these matrices are linearized 
about the steady deflection, the inertia and stiffness matrices are symmetric whereas 
the damping matrix is antisymmetric; but the nonconservative aerodynamic terms 
will make the damping and stiffness matrices asymmetric. 

Inertia Matrix 

The element inertia matrix for the ith element is written as 


where 


[A/(l 


[A/uv] (A/vw) [Afy^] 
[Mftiv] (A/uiw| [A/u;^] 

.(Afi^] [A/^u,] 


{C.l) 





[Mvu>] “ (01 



meg sin dx^ 


JO 
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mefCO»io{H){H)^dx{ 


[Mn\- l‘' mkl{H){Hfdxi 


Damping Matrix 

The element damping matrix for the ith element is written as 



•(Owl 


[^V^l 

[<^.•1 - 

[C/wv] 

(Cu^u,) 





[<^hV 


wheTO 



2megcos9o 


dXf 



-2m^p{H){H}^ - 2mtgs\n9o - {H){Htf 


dzi 


[OvA = lo| 


(Cu;u,] =* [0] 
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W - 10) 


M - 101 


Sf,i|Tness Matrix 


m 


where 


. l^^l [^'•^1 


(C.3) 




lK^\-f'{l{Eh-El,-EA.t),intlo{H"}{H’'f 

® ^ ji n 
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[K^^] - * j^(jS7^ea^A - ^ 82)60 cos 

+ fo<a sin^o{^^^}{J^^}^ 
-E(72sm»o{H"}|j^fj} 

+ tnuj sin iD[H] 

— m«2asintfo{^^^){^^}^ 

- (E/, - B/,) sin 2»o{H") {E")’’(<>o) {ff# )’’ 

+ {El, - E/,)cos2>o{H"){H"}’'{”’o){«*)'' 

+ oy{E"){E'}’'{«5){wV}1<<*.- 


[E„^] — ' f(BAsa*^ - EE2)«'osin«o{ff")|Hi} 

— Fo«a COS 0o{H’^} {8[,f>}^ 

T* 

+ EC?2Cos«o{H')|wJ| 

J1 

+ megXCOs$o[H^}{H^} 

+ (Elg - EIy)Q092$(i{H'^}{H'*}^{vQ}{H^}^ 
-F {Elg - EIy)B\n20o{H'^}{H'*]^{wQ}{H^}^ 

+ G/{E'){E*}’'{t-})|EV} H*.- 


\K^\ - jf ' HG J + Eo*^ - EAk\4 + EBi 

+ec,|e5||e;}^ 

+ - k%^yos28o{H^}{H^} 




dxi 


The column vectors {H}, {H^}, and {vg} are defined as 
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{«)’■- IHl Hi Hi H,i 


Wf-IH^ Htii 


*■ l*'Ol **01 *^02 *'02! 


A similar definition holds good for (»}) ax** 



Also ^0 is defined as 


#0 “ ^ + ^0 
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APPENDIX D. EQUIVALENT-BEAM PROPERTIES 

Mathematical expressions of equivalent properties for modelling twin beams 
(Fig. 17) as an equivalent beam are presented here. 


(EA), = {EA)i + (£7A)2 

(El,), « (£?/.), + {EI,)i + ViiEA), + tiUEA)2 

(^^)2 

me * mi + m2 


The subscripts and 2 refer, respectively, to the equivalent beam, flexbeam one, 
and flexbeam two; and Ly is the length of the flexure. These expressions become 
simpler for the present case where the two beams are identical. 
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Table 1. Values of parameters of imiform blade for numerical results 


Parameter 

Numerical Value 


0.014480 1- 1.15) 

ErMimd^n* 

0.026655 (w, - OJ) 
0.166008 (h>» - 1.5) 


O.OQ0025(o;^ » 2.5) 
0.005661(0^^ - 5) 

*»../« 

0 

kmjR 

0.025 


1.5 

ejR 

ir/40 

Pf 

0.05 rad 

or 

0.1 

7 

5 

a 

6 


0 

Cl 

6 

do 

0.0005 

dx 

0 

<h 

0 


0 

kk 

1.15 










Table 2. Steady tip deflections of a biogeless rotor blade 


Number of 
elemeote 

Vo, JR 

vfo,jR 


2 

•0.00350 

0.00487 

•0.0438« 

3 

•0.00342 

0.00459 

•0.04334 

4 

•0.00338 

0.00446 

•0.04315 

5 

•0.00337 

0.00440 

•0.04307 

0 

•0.00330 

0.iIK)436 

-0.04302 

7 

•0,00335 

0.00434 

•0.04299 

8 

•0,00335 

0.00433 

-0.04297 


(Cr/<r ■■ 0.1, ■■ 0.05 rad, 7 ■■ 6.0, w, — 1,50, Uw 1.150, ■■ 2.50) 
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Tabl^ S. Fundamental Coupled natural frequencies of a hingeless rotor blade 


Number of 
elemeots 

(Jv/n 

uTw/n 

vtia 

f% 

m. 

1.6100 

1.1251 

2.5105 

3 

1.5184 

1,1221 

2,4872 

4 

1.5182 

1.1214 

2.4787 

5 

1.5181 

1.1212 

2.4746 

0 

1.5181 

1.1211 

2.4724 

7 

1.5181 

1.1211 

2.4710 

8 

1.5180 

1.1210 

2.4702 


(Crier ■■ 0,1, fip 0.05 rad, 7 ■■ 6.0, ■■ i.Sn, (tjy, ■■ l.lSfl, ■■ 2.50) 
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Table 4. Real pp;rts of the flutter eigenvalues of a hingeless rotor blade 


Number of 
modes 

Lead'lag 

Flap 

ToriioD 

3 

•0.03074 

-0.31488 

•0.35148 

4 

-0,03040 

-0.31383 

•0.35040 

5 

•0.03034 

•0.31443 

•0.35207 

0 

•0.03034 

•0,31442 

>0.3520e 

7 

-0.03034 

-0,31448 

-0.35200 


(Cr/fl^ " 0.1, "" 0.05 rad, 7 ■■ 6.0, u;« ■■ 1.50, w» — 1.150, ■■ 2.50) 
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Table 5. Steady tip deflections of an articulated rotor blade 


Number of 
elements 

Vo, JR 

V’OujR 

A 

2 

.0.06040 

0.01274 

.0.04476 

3 

.0.06001 

0.01238 

.0.04425 

4 

4).05987 

0.01224 

.0.04405 

5 

.0.05080 

0.01217 

.0.04396 

6 

4).05977 

0.01212 

.0.04391 

7 

.0.05975 

0.01210 

.0.04388 

8 

4).05973 

0.01208 

.0.04386 


■■ 0.1, ■« 0.05 rad, 7 — 5.0, Hinge offset ■■ 0.06A, “ 2.50) 




- 1 - 














Table 6. Fundamental Coupled natural frequencies of an articulated rotor blade 


Number of 
elements 


k>w/n 


2 

0.2909 

1.0440 

2.5285 

3 


1.0440 

2.5049 

4 

0.2099 

1.0440 

2.4963 

5 

0,2900 

1.0440 

2.4923 

6 

0.2999 

1.0440 

2.4001 

9 

7 

0.2990 

1,0440 

2.4889 

8 

0.2999 

1.0440 

2.4878 


(Cr/<r ^ 0.1, ffp ■■ 0.05 rad, 7 ■■ 5.0, Hinge offset ■■ 0.06/7, UI 4 -■ 2.5fi) 
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Table 7. Real parts of the flutter eigenvalues of an articulated rotor blade 


Number of 
modes 

Lead-lag 

Flap 

Torsion 

3 

-0.00055 

-0.33926 

-0.39160 

4 

-0.00051 

•0.34108 

•0,39338 

5 

-0.00957 

-0.34498 

-0.39471 

6 

-0.00953 

•0.34232 

-0.39448 

7 

•0.00953 

-0.34230 

-0.39449 


{Crfo ■■ O.l, /Jj ■■ 0.05 rad, 'f »■ 5.0, Hinge offset ■■ 0.06H, » 2,50) 










Table 8. Structural properties of a bearingless blade 


Property 

Outboard blade 

Each flexbeam 

Torque tube 


0.026655 

0.013328 

0.001 

0.014486 

0.007243 

0.001 

EIufmon^R^ 

0.166908 

0.083454 

0.001 

GJlmon^R* 

0.000925 

0.000463 

0.02 


0.005661 

0.002830 

0.02 

kmijR 

0 

0 

0 

km^lR 

0.025 

0.025 

0.01 

^Al^m 

1.5 

1.5 

1 

mimo 

1 

0.5 

0.1 
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Table 0. Steady tip deflections of a twinbeam model 
(Zero inboard pitch) 


Number 

Flexbeam 

one 

of 

Flexbeam 

two 

elements 

Outboard 

blade 





1 

1 

2 

0.001613 

^.003750 

0.011214 

•0.030734 

2 

2 

1 

0.001614 

•0.003704 

0.011280 

•0.030634 

2 

2 

2 

0.001613 

•0.003750 

0.011213 

•0.030735 

2 

2 

3 

0.001613 

•0.003747 

0.011206 

•0.030732 

3 

3 

— j 

2 

0.001613 

•0.003750 

0.011213 

-0.030735 


(C*r/<r "" 0.1, 0.05 rad, 7 ■■ 5.0, u>^ ■■ 1.87(1, ■■ 1.15(1, u>^ ■■ 2.01(1) 














Table 10. Fundamental natural frequencies of a twinbeam model 

(Zero inboard pitch) 


Number 

of 

elements 




Flexbeam 

one 

Flexbeam 

two 

Outboard 

blade 

a;«/0 


w^/0 

1 

1 

2 

1.77073 

1.15018 

2.88149 

2 

2 

1 

1.78825 

1.15027 

2,89491 

2 

2 

2 

1.77962 

1.14977 

2,88144 

2 

2 

3 

1.77877 

1.14971 

2,88094 

3 

3 

2 

1.77961 

1.14974 

2.8S141 


(Cr/o 0.1, fip 0.05 rad, 7 5.0, » 1.87(1, » 1.150, >■ 2.010) 












Table 11. Real parts of the flutter eigenvalues of a twinbeam model 

(Zero inboard pitch) 


Number of 
modes 

Lead'lag 

Flap 

Torsion 

3 

0.00506 

•0.35029 

-0.38113 

4 

0.00450 

-0.35228 

•0.38094 

5 

0.00445 

-0.35267 

-0.38119 

0 

0.00445 

-0.35267 

-0.38119 


(Cr/<T — 0.1, — 0.05 rad, 7 — 5.0, u, — 1.87(1, — 1.15(1, U 4 , » 2.91(1) 












Table 12. Fundamental natural frequencies at various thrust levels 

(Zero inboard pitch) 








LAG HINGE 



Fig. 1 Schematic of an articulated rotor hub and root. 



BLACK AMD WHITE PHOTOGRAEK 



Fig. 2 Bearingless Main Rotor (BMR) 
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Fig. 3 Flexure of BMR 
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AXIS OF 
ROTATION 



Fig. 4 Blade coordinate systems and deflections 
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Fig. 5 Deformed-blade-section aerodnamics 





Fig. 6 A finite element showing nodal degrees of freedom 

. 1 
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Fig. 7 (a) Flap bending (b) lead-lag bending (c) torsion. Steady 
tip deOections of a stiffdnplane hingeless rotor blade; = 
l.sn, Wu, *= l.isn, = 5 . 0 n, 7 = s o, = 0.05 rad. 
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I 

I ■ 

Fig 8 (a) Flap bending (b) lead-lag bending (c) torsion. Steady i 

tip deflections of a soft-inplane hingeless rotor blade; , 

0.70, *“ I'isn, =* 5.00, If ■= 5.0, 0.05 rad. i 
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Fig. 10 Stability boundaries for a hingeless rotor blade; Uui 
1.150, — 2.50, 7 — 5,0, Pp — 0.05 rad. 
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PRESENT 

EXPERIMENT 



Fig. 1 1 Comparison of hingelesss rotor blade results with experimen 
tat results. 


f; 

i 




Fig. 12 a) Flap bending, b) lead-lag bending c) torsion. Steady 

tip deflections of an articulated rotor blade; hinge oiTset=0.06R, 
» 2.sn, 7 « 5.0, — 0.05 rad. 




ORfGiiVAsl 




POOR QUAU1Y 


2.0 


2.4 


2.2 


2.0 


1.8 


1.0 


1.4 


E 


1.2 


1.0 


.6 


.4 


.2 


I TORSION 


SION^ 


C^o 

I > 

^ .1 Jt ,3 


FLAP 




J ^ 


-.0 -.5 -.4 


*.3 
RC \ 


LAG 


-.2 -.1 


Fig. 13 Root locus plot for an articulated rotor blade; hinge offset= 0 . 06 R, 
hj^ 2.50, 7 « 5.0, * 0.05 rad. 
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Fig. 14 Analytical model of a bearingless blade. 



Fig. 15 Refined finite element for multiple-load-path blades. 
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Fig. 16 Pitch distribution of a bearingless rotor blade. 




Fig. 17 Twinbeam multiple-load-path blade model. 
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Fig, 18 Root locus plot of lead-lag mode of stiff-inplane, soft-in- 
torsion twinbeam model (zero inboard pitch); Wv -» 1.870, 
o^u, » l,isn, ■■ 2,010, ^ ■■ 5,0, ■■ 0,05 rad. 
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Fig. 10 Root locus plot of lead>lag mode of stifT-inplane, soft- 
in*torsion twinbeam model (zero inboard pitch; coupling 
structural stiffness terms of nonlinear origin suppressed 
in flutter analysis); — 1.870, uw — U5f], — 2.910, 

7 » 5.0, *a C,u3rod. 
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Fig. 20 Root locus plot of lead*lag mode of soft-inplane, soft-in- 
torsion twinbeam model (zero inboard pitch); utv “ 0,900, 
Uw - 1.150, - 2.910, If - 5.0, - 0.05 rad. 
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Fig. 21 Root locus plot of lead-lag mode of soft-inplane, soft- 
in-torsioD t^vinbeam model (zero inboard pitch; coupling 
structural stiffness terms of nonlinear origin suppressed 
in flutter analysis); wv ™ 0.900, uw ■* 1.150, »» 2.910, 

7 5.0, Pp »» 0.05 rad. 
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Fig. 22 Root locus plot of lead-lag mode of stiff-inplane, stiff-in- 
torsion twinbeam model (zero inboard pitch); a;„ = 1,870, 
“ l.lSfl, 0 ?^ ■* 5.380, 7 5.0, « 0.05 rod. 
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Fig. 23 Root locus plot of lead-lag mode of stiff-in plane, stiff-in- 
torsion twinbcam model (zero inboard pitch); coupling 
structural stiffness terms of nonlinear origin suppressed 
in flutter analysis); — 1.870, Wu, ■■ 1.150, — 5.380, 

-7 5.0, « 0.05 rad. 
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Fie 24 Root locus plot of lead-lag mode of f J 
in-torsion twinbeam model (same pitch for inboard and 
outboard); wv - 1-870, a^u. - 1-150, - 2-910, 'i - 5,0, 

— 0.05 rad. 
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Fig. 25 Root locus plot of lead-lag mode of soft^-inplane, soft- 
in-torsion twinbeam model (same pitch for inboard and 
outboard); » 0.90n, — 1,150, — 2.910, -y — 5.0, 

/^p— 0.05 rad. 
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Fig. 26 Root locus plot of lead-lag mode of stiff-inplane, stiff- 
in-torsion twinbeam model (same pitch for inboard and 
outboard); ~ 1.870, Uui -■ 1.150, 5.380, 7 5.0, 

ftf ■« 0.05 rad. 




